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Abstract 



o 

■ The fully relativistic theory of the Zeeman splitting of the (ls) 2 2s hyperfine-structure 

levels in lithiumlike ions with Z = 6 — 32 is considered for the magnetic field magnitude in 
the range from 1 to 10 T. The second-order corrections to the Breit - Rabi formula are cal- 
Q ■ culated and discussed including the one-electron contributions as well as the interelectronic- 

inter action effects of order 1/Z. The 1/Z corrections are evaluated within a rigorous QED 
| approach. These corrections are combined with other interelectronic-interaction, QED, nu- 

clear recoil, and nuclear size corrections to obtain high-precision theoretical values for the 
Zeeman splitting in Li-like ions with nonzero nuclear spin. The results can be used for a 
precise determination of nuclear magnetic moments from g-factor experiments. 



q ; PACS number(s): 32.60.+i, 31.30.Gs, 31.30.Jv, 12.20.Ds 

o 1 Introduction 



CZ2 

J^'. Recent measurements of the g factor of hydrogenlike carbon and oxygen have reached an accuracy 
o I of about 2 • 10~ 9 [1-3]. The experiments were performed on a single hydrogenlike ion confined 
in a Penning trap with a strong magnetic field (B = 3.8 T). These measurements considerably 
stimulated theoretical investigations of this effect [4-22] . Besides a new possibility for tests of the 
magnetic sector of quantum electrodynamics (QED), these investigations have already provided 
^ \ a new determination of the electron mass (see Refs. [3,23] and references therein). Extensions 
^ \ of these experiments to systems with higher nuclear charge number Z and to ions with nonzero 
nuclear spin would also provide the basis for new determinations of the fine-structure constant 



[8,24,25], the nuclear magnetic moments [24], and the nuclear charge radii. 

Extending theoretical description from an H-like to a Li-like ion, one encounters a serious 
complication due to the presence of additional electrons. A number of relativistic calculations of 
[ the g factor of Li-like ions were carried out previously [26-30]. However, to reach the accuracy 
^ ' comparable to the one for H-like ions, a systematic quantum electrodynamic (QED) treatment is 
^ ■ required [16,31-35]. 

For both H- and Li-like heavy ions with nonzero nuclear spin the ground-state Zeeman splitting 
caused by the magnetic field in the range from 1 to 10 T is much smaller than the hyperfine split- 
ting and, therefore, the consideration can be conveniently reduced to the g factor value [19,35]. 
However, for H-like ions with Z = 1 — 20, which are being under current experimental investiga- 
tions at Mainz University, the Zeeman splitting is comparable with the hyperfine splitting if the 
magnitude of the homogeneous magnetic field does not exceed 10 T. This requires constructing 
the perturbation theory for degenerate states. To a good accuracy, the solution of the problem is 
given by the well-known Breit - Rabi formula [36-39]. The aforesaid experimental precision has, 
however, shown the need for an improvement of the Breit - Rabi formula for H-like ions [21]. 
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In the present paper, we consider the Breit - Rabi formula for the 2s hyperfme-structure levels 
in lithiumlike ions. Evaluations of the coefficients of this formula should include corrections de- 
pending on the nuclear g factor. Besides a simple one-electron lowest-order nuclear-spin-dependent 
contribution, one should also calculate the second-order corrections caused by the hyperfine inter- 
action and the interaction with the external magnetic field, taking into account the presence of 
the closed (Is) 2 electron shell. We perform such calculations in the range Z = 6 — 32, where the 2s 
HFS splitting can be comparable with the Zeeman splitting if the magnitude of the homogeneous 
magnetic field is in the range under consideration. The calculations are based on perturbation 
theory in the parameter 1/Z within a rigorous QED approach. The contributions of zeroth and 
first orders in 1/Z are taken into account for the magnetic-dipole correction and the contribution 
of zeroth order is taken into consideration for the electric-quadrupole correction. Also, the in- 
dependent correction is calculated, including the contributions of zeroth and first orders in 1/Z. 
The obtained results are combined with other corrections to get accurate theoretical predictions 
for the Breit - Rabi formula coefficients for lithiumlike ions with nonzero nuclear spin. These pre- 
dictions will be important for experimental investigations that are anticipated in the near future 
at University of Mainz and GSI [40]. 

Relativistic units (h = c = 1) and the Heaviside charge unit (a = e 2 /4n, e < 0) are used in 
the paper. In some important cases, the final formulas contain h and c explicitly to be applicable 
for arbitrary system of units. 



2 The Breit — Rabi formula in the lowest-order one-electron 
approximation 

We consider a lithiumlike ion with nonzero nuclear spin J in a state of the valent electron with 
the total electron angular momentum j = 1/2. The ion is placed in a homogeneous magnetic 
field B directed along the z axis. The magnetic splitting is linear with respect to B only if 
one of the following conditions is fulfilled: either AE mag <C AEhfs or AE mag ^> AEhfs, where 
AEhfs = E(F + 1) - E(F), E{F) = E nK + e Ms (F), F = I ± 1/2 is the total atomic angular 
momentum, and Ehfs(F) is the hyperfme-structure shift of the valent electron Dirac state with the 
one-electron energy in the Coulomb field of the nucleus 

En* = N m e . (1) 

Here n is the principal quantum number, k = (— iy +l+ ^{j + |), I — j ± \ defines the parity of 
the state, n r = n — \k\ is the radial quantum number, 7 = k 2 — (aZ) 2 , N = \Jn 2 . + 2n r 7 + k 2 , 
and m e is the electron mass. It should be emphasized that in case the second inequality is 
fulfilled AE mag must be much less than the distance to other Dirac's levels. In the intermediate 

— * 

B case, AE mag ~ AE'hfs? we must take into account mixing the HFS sublevels with the same 
Mp, where Mp = —F, —F + 1, F — 1, F is the z projection of the total angular momentum. 
For the states with the total electron angular momentum j = 1/2, there are only two HFS levels 
F = 1-1/2 and F' = 1+1/2 with the same M F = -7+1/2, ...,1-1/2. This greatly simplifies the 
theory. In what follows, we restrict our consideration to the ground state of the valent electron. 
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Denoting AE mSLg = E — — — one can derive f or the Zeeman splitting 



AE ma Jx) = AE- 



(2s) 
HFS 



AMp 

aiM F x ± -\ l-\ 

2 V 27+1 



where re = fj, Q B/AE^p S , /i = |e|ft/(2m e c) is the Bohr magneton, 

ai = -g'i , 
ci = Qj + , 

c 2 = + g'i) 2 



(2) 



is the ground-state bound-electron g factor of the lithiumlike ion, 

9j = 9d + Ag hlt + As-qed + + A^ns + A^np 
<7d is the one-electron Dirac value for a point-charge nucleus, 

2[ V /2T2^+1] (aZ) 2 
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+ ..., 



(3) 
(4) 
(5) 



(6) 



(7) 



7 = a/1 — (aZ) 2 , Ag int is the interelectronic-interaction correction, Agq ED is the QED correction, 
is the nuclear recoil correction to the bound-electron g factor, A^ns is the nuclear size 
correction, A^np is the nuclear polarization correction, g\ is the nuclear g factor expressed in the 
Bohr magnetons, 

— foz + AyW), (8) 



g'i 



is the proton mass, gi = /i/(/iyv/) , ji = (II\fi z \II) is the nuclear magnetic moment, \i z is 
the z projection of the nuclear magnetic moment operator jl acting in the space of nuclear wave 
functions \IMj) with the total angular momentum / and its projection Mj, fj,^ = \e\h/(2m p c) 
is the nuclear magneton, and Ag^cc is the recoil correction to the bound-nucleus g factor (see 
section[5]). Eq. (j2|) is usually called the Breit - Rabi formula (see, e.g., Refs. [21,36-39]). It covers 
Zeeman splitting from weak {x <C 1) to strong (x ^> 1) fields including the intermediate region. 
For F' = I + ~ and Mp = ±(J + |) the splitting is linear in the first order of perturbation theory 
under arbitrary magnetic induction, 



A£ mag (x) = AE 



(2s) 
HFS 



- ± dix 



where 



di = \gj - ig'i 



(9) 



(10) 



1 In the present paper, the energy of a Zeeman sublevel AE mag is counted with respect to the mean energy 

— - of the hypcrfinc structure doublet [38,39]. To count the energy from the hyperfme centroid of 

the doublet [36,37], one should use the relation 



AE, 



AE- 



(2s) 
HFS 



mag 



2(27+1) 
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and the "— " and "+" signs refer to Mp = —(I + |) and Mp — I + ~, respectively. 

For Li- like ions with I = 1/2, F = and F' — 1 and, therefore, the two mixed sublevels have 
Mp = 0. In this case the Breit - Rabi formula takes the form 

AF mag (x) = ±^ V / TW, (11) 

and for M F = ±1 the effect is described by Eq. (J5J) with d\ = \{gj ~ <?/)• 

If the magnetic field is strong, AF mag ^> AE^p S , Eqs. (EJ), (EJ), and fTTTT) convert into formulas 

for the anomalous Zeeman effect of the 2s state. 

In case the energy-level shift (splitting) due to interaction with B is much smaller than the 

hyperfine-structure splitting, AF mag <C AE^p S , we can express the linear-dependent part of this 

shift in terms of the atomic g factor, 

AE (2s) 

AF mag = + g{F)^BM F , (12) 

where, to the lowest-order approximation (see, e.g., Ref. [38]), 

77? 

g(F) = g D Y cl (F) - -^ 9l Y^ c (F) , (13) 



Y(F\- F{F + 1) + 3/4 - 1(1 + 1) _ J - ^ forF = /-I 
Cl[) ~ 2F(F + 1) -\^- forF = / + i' 1 ' 

Y U) (F) _ F(F + 1) + J(J+ 1) - 3/4 _ forF=J-| 
nucl J " 2F(F + 1) forF = J+|' 1 J 

The total one-electron 2s ^-factor value of a Li-like ion with nonzero nuclear spin can be 
represented by 

g(F) = (g D + Ag int + Ag QED + Ag£ c + Ag NS + Ag NP )Y cl (F) 

-^( 9I + Ag^)Y^(F) + 8g$ s (F) , (16) 

where the HFS correction 5g^p S (F) = Sg^p S ^(F) + ^<7hfs(Q)(-^) 1^5] is briefly discussed below. 

3 Hyperfine-interaction corrections to the ground state g 
factor 

Let us start our consideration of the HFS correction to the ground-state g factor of a Li-like ion 
with the one-electron approximation. In this approximation, the interaction of the ion with the 
magnetic field can be represented as 
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Here vl describes the interaction of the valent 2s electron with the homogeneous magnetic field, 

V^ = -e(a-A)= l -^(a-[Bxr\), (18) 

where the vector a incorporates the Dirac a matrices, and 

V$* = -{jil.$) (19) 

describes the interaction of the nuclear magnetic moment \x with B. The hyperfine-interaction 
operator is given by the sum 

VhFS = V^fs + ^HFS ) ( 20 ) 

where V^pg and V^p S are the magnetic-dipole and electric-quadrupole hyperfine-interaction oper- 



ators, respectively In the point-dipole approximation, 



HFS ~ 4tt r 3 ' { ' 



and, in the point-quadrupole approximation, 



m=2 



^h?s = -« E Q^m(n). (22) 

m=-2 

Here Q 2 m = X)^=i ffC^m^i) is the operator of the electric-quadrupole moment of the nucleus, 
V2m — C2 m (n)/r 3 is an operator that acts on electron variables, ft — f/r, fij = f^/rj, r is the 
position vector of the electron, is the position vector of the i-th proton in the nucleus, Ci m = 
a/47t/ (21 + 1) Yi m , and Yi m is a spherical harmonic. It must be stressed that the electric-quadrupole 
interaction should be taken into account only for ions with I > 1/2. 

An unperturbed atomic eigenstate that corresponds to given values of F and Mp is a linear 
combination of products of electron and nuclear wave functions, 

\nl 3 IFM F )= C^lnljm^IMr). (23) 

rrij ,Mt 

Here C^J Mi are the Clebsch-Gordan coefficients, \nljrrij) are the unperturbed one-electron wave 
functions, which are four-component eigenvectors of the Dirac equation for the Coulomb field, 
with the total angular momentum j and its projection rrij. 

In the one-electron approximation, the magnetic-dipole and electric-quadrupole hyperfine- 
interaction corrections to the ground-state g factor of the Li-like ion are given by 

r onc-el.(2s) _ 2 ^ ST^ r FM F r FM F ,j M \ B l"A ra I^HFS \ U I , 

rrijMim'.M'j n 

(24) 

where \v) = \20^rrij) and \v') = |20|m^) are the 2s states of the valent electron with the angular 
momentum projections rrij and rrij, respectively, |n) = \nljrrij), e v = E2-1, and e n = E nK . The 
summation in f[2~4l) runs over discrete as well as continuum states. The corresponding diagrams 
are presented in Fig. [H 
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The total hyperfine-interaction correction to the ground-state g factor of the Li-like ion is given 

by 

S 9%f S = ^HFS(m) + ^HFS(Q) ( 25 ) 

with 

^HFSM = a2Z h^ p l Y ^ {F)[S2{aZ) + + h C ^ Z) + ' ' ' ] (26) 

and 

^hfs(q) = a ' z3 ^Q (^) Vw>(F)[r a (oZ) + ±B Q {aZ) + j~ 2 C Q (aZ) + ...]. (27) 
Here the angular factor is 

F n («)(F) = J ^-D(2m) Ior ^ 2 (28) 

I 27TT forF = / + i 

and Q = 2(II\Q2q\II) is the electric-quadrupole moment of the nucleus. The functions 

s ^> %^;W) ^ r (29) 



and 

T ^ az ) = - ~,~is2„(o^,^ ( 3 °) 



2160 one-el. (2s) 

23a*Z3Q(^) 2 Y£\F) ^ HFS(Q) 



determine the one-electron contributions, which are discussed in detail in Ref. [19]. For the point- 
charge nucleus, the functions ^(aZ) and T 2 (aZ) are [19,35] 



S 2 (aZ) 



±1. 



N + 



io(iv + r 

3~/V 



+ 



{aZf f 2(iV+ 1) 
7(7 + 1) [3-4:(aZy 



+ 1 



3N { N + 2 

l + ^(aZ) 2 + ... (31) 



144 
and 

192[(iV + 7 + 1){18 + 24 7 - 12iV + 87iV 2 } + 15(1 + 7)] 



T 2 (aZ) 



237iV 3 [15 - 16(aZ) 2 ](iV + 7 + l) 2 

427 

1 + 276 (aZ)2 + '--' (32) 



where N = y/2(l + j). 

The interelectronic-interaction correction B^(aZ) can be calculated within the rigorous QED 
approach [35]. The interaction of the electrons with the Coulomb field of the nucleus is included 
in the unperturbed Hamiltonian, i.e. the Furry picture is used. The perturbation theory is 
formulated with the technique of the two-time Green function (TTGF) [41,42]. To simplify the 
calculations, the closed (Is) 2 shell is regarded as belonging to a redefined vacuum. With this 
vacuum, the Fourier transform of TTGF can be introduced by 



00 



Q{E- x'\ x)5(E -E') = — [ dx°dx'° exp(iE'x'° - iEx°) 

2m J 



(33) 

-CO 

x (0 {u) 2\T^(x' ,x')^(x°,x)\0 {ls) 2), 
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where ip(x°, x) is the electron-positron field operator in the Heisenberg representation and T is 
the time-ordered product operator. The energy shift of a state a can be expressed in terms of the 
TTGF defined by 

g aa {E) = {u a \Q(E)\u a ) = J dxdx'ul(x')G(E;x';x)u a (x), (34) 

where u a {x) is the unperturbed Dirac wave function of the state a. Using the Sz.-Nagy and Kato 
technique [43], one can derive for the total energy shift AE a = E a - E [ a ] [41,42] 

£-JdEAEA 9(m (E) 

AEa= l + r £-JdEAgUE)> (35) 
r 

where AE = E-E$\ Ag aa (E) = g aa (E)-g^ (E), and (£) = (E-E^)' 1 - The integrals in the 
complex i?-plane are taken along the contour V which surrounds the pole of g aa (E) corresponding 
to the level a and keeps outside all other singularities. The contour T is oriented counter-clockwise. 
To first three orders of the perturbation theory, the energy shift is given by 

AE °' = ^Ti£ dE AEA,J ™ {E) ' (36) 



ae «' =hl iEAEA ^ E) - {^l iEAEA ^) {■Li iEA ^ E) 



(37) 



+ {hl dEAEA ^w) (hl iEA ^ E) ) 2 ■ < 38) 

The redefinition of the vacuum changes iO to — iO in the electron propagator denominators cor- 
responding to the closed (Is) 2 shell. In other words it means replacing the standard Feynman 
contour of integration over the electron energy C with a new contour C (Fig. [2]). The second- 
order contribution is defined by the diagrams presented in Fig. [TJ Its evaluation according to Eq. 
(13*7)) yields formula (T241) . In the formalism under consideration, the lowest-order interelectronic- 
interaction and the radiative corrections to Eq. ()24l) are described by the third-order diagrams 
presented in Fig. [3] and, according to Eq. ( 1381) . by some products of the low-order diagrams 
depicted in Figs. H] and [5j According to Fig. El to separate the interelectronic-interaction correc- 
tions, the contour C must be divided into two parts, C and C int . The integral along the standard 
Feynman contour C gives the one-electron radiative correction. The integral along the contour 
Ci n t describes the interaction of the valent electron with the closed shell electrons. Formula (13*51) 
allows one to evaluate the interelectronic-interaction correction B^(aZ) [35]. The results of this 
evaluation will be presented in the next section together with other related corrections to the 
Breit - Rabi formula. 
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4 Corrections to the Breit — Rabi formula for the ground 
state 

Now we assume that the Zeeman splitting AP mag of the 2s HFS levels F = I — 1/2 and F' = 
I + 1/2 is much smaller than the distance to other levels but is comparable with AE^ S . The 
unperturbed eigenstates form a two-dimensional subspace f2 = {ll^ -*), |2^°- ) )}, where = 
\20±IFM F ), |2(°>) = |20| IF'M F ). Employing the perturbation theory for degenerate states [42] 
with energy e v we denote the projector on Vt by 

p(°) = ^(°))( z -(°)|. (39) 

i=i 

We project the Green function Q(E) on the subspace Q 

g(E) = pV>g(E)P® , (40) 

where, as in Eq. (1341) . the integration over the electron coordinates is implicit. In this case we 
can choose a contour T in the complex P-plane in a way that it surrounds all g(E) poles, which 
correspond to the states under consideration, and keeps outside all other singularities of g(E). As 
in the case of a single level, to the zeroth-order approximation one easily finds 

We introduce the operators K and P by 

K = ^-jdEEg(E), (42) 



P =^l% dEa[E) - <43) 
As it is shown in Ref. [42], the energy levels are determined from the equation 

det(#- E) = 0, (44) 

where 

H = P^KP-^ . (45) 
The operators K and P are constructed by formulas (142 j) and (143jl 

K =K<® + + i^ 2 ) + K® + . . . , (46) 

p =P (0) + p(D + p(2) + p(3) + _ _ ? (47) 

where the superscript indicates the order of the perturbation theory in a small parameter. The 
operator H is 

H = PT (0) + + # (2) + # (3) + . . . , (48) 



8 



where 

if <°> = , (49) 
H W = K (D _ i P (i)#(o) _ 1^(0) p(D ; (50) 

H (2) = K (2) _ I p(2) K (0) _ l K (0) p(2) _ 1 p(l) A Kl) _ p(l) 

2 2 2 2 

+ ^ pW p(l)^(0) + 3^(0) p(l) p(l) + 1 p(l) K (0) p(l) ^ / 51 N 

8 8 4 

p/(3) = ^(3) _ 1 p(3) j^(0) _ 1^(0) p(3) _ 1 p(l) j^(2) _ 1^(2) p(l) _ 1 p(2) ^(1) _ 1 p(2) 

2 2 2 2 2 2 

+ - pi 1 ) p(2)^(0) , 3^(0) p(l) p(2) , 3 p(2) p(l) K (0) , 3^(0) p(2) p(l) 

8 8 8 8 

+ - P^K {0) P (2) + - P^K {0) P (1) + - P (1) P^K {1) + -K (1) P (1) P (1) + - P (1) if (1) P (1) 
4 4 8 8 4 

_ lp(i) p(i) p(i)^(o) _ AirW pW P« P« 
16 16 

_ A p(i) p(i) 7^(o) p(i) _ A p(i) &:(o) p(i) p(i) f ko 1 ) 

16 16 V ; 

Taking into account only the relevant contributions of kind a x fi/fi^ x -B and a x B x B, 
where a comes from the interelectronic interaction, we obtain for the third-order term in Eq 

= ds 1 d£ A£A9 * < £ > 4 1 [(i i d£ As "' < £ i d£ A£A9 - 21 < £ >) 
+ (ii d£A£A ^( £ ») (^// ba ^ 



(53) 



where AE = E — e v , j, k = 1,2. 

Keeping only the three lowest-order terms in B, we get the following equation for the perturbed 
energies: 



ho(F) + h 1 {F)B + h 2 {F)B 2 - E h^F, F')B + h 2 {F, F')B 2 
hiF', F)B + h 2 (F', F)B 2 ho(F') + h^B + h 2 (F')B 2 - E 



Here F = I - ±, F' = I + \, 



is the energy of the HFS level, 



. (54) 

h (k) = E(k) (55) 



hi(k) = ^E ( ^(k, k) + AE$(k, k) + AE®{k, k) + AE$(k, k) + AE®(k, k)] 

+ {Ag int + Ag QED + AgQ + Ag m + Ag^ P )Y el (k)fi Q M F - Ag^Y^ {k)fi N M F (56) 
= g(k)fM)M F , 
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h 2 (k) = -^[^(k, k) + AE^(k, k)} , (57) 

+ (A int + A QED + Ag> + A NS + A NP )// - A^Vjv , (58) 

fcfr fc) = ^ fc) + AE®(j, fc)] , (59) 

where j, fc = F, F'. A int , Aqed, Arec, Ans, and Anp are the interelectronic-interaction, QED, 
nuclear recoil, nuclear size, and nuclear polarization corrections. They are similar to the corre- 
sponding corrections to hi(k) but have a different angular factor as well as A^i It should be 
noted that we have neglected here terms describing virtual transitions into excited nuclear states 
via the direct interaction of the nucleus with the magnetic field [22]. The energy shifts are 

A£ §)0'> fc )=EE ^^m;(^i^W) > ( 60 ) 

mi Mj m'.M' 

J 3 i 

AE m -Aj,k)= Y V Ci^^C^FJlMAlW MM'), (61) 

£ o>^, m c^ • (62) 

where 

% 2fZ Ufl '_ ' (63) 
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T (3a) 
fj,,Q,B 



f E ( E 



\E V £ni)\£v ^712) 



(vlV^lnJimlW^inzcllWlv'c) 



+ (vlV^ln^MlWlnz^inilWlv') + ^|iy|n 1 )(n 1 |V^ e) |n 2 )^2c|/(0)|^c) 



E E 



(v\vi e) \n)(n\W\v)(vc\I{0)\ v'c) 



+ (v\Vi e) \n)(nc\I(0)\vc)(v\W\v') + (v\Vi e) \v)(v\W\n)(nc\I(0)\v'c) 



(64) 



r(36) 



E 

711,712 



^|V^ e Vi)KI^I^)(ri 2 c|/M|ct/) 



(e„ - e ni )(e w - £n 2 ) 
+ (^|Vl e) |ni)(«ic|7(a;)|cn 2 )(n 2 |W|^} + M^|r2i)(r2 1 |yl e V2)(n 2 c|/(^)|c^) 



E E 



(v\vi e) \n)(n\W\ v) (vc\I (cu) \ cv') 



+ (v\vi e) \n)(nc\I(u)\cv)(v\W\v') + (v\vi e) \v) (v\W\n) (nc\I(u)\cv') 



E E 



(v\vi e) \n)(n\W\v)(vc\l'(u;)\cv') 
+ {v\vi e) \n)(nc\I\uj)\cv}(v\W\v'} + (v\Vi e) \v)(v\W\n)(nc\I'(oj)\cv') 

B B 

+ J2 ^^\V^\v)(vc\I'\u)\cv)(v\W\v') 



(65) 



T (3c) _ 



f E ( E 



{e v £ ni ){£ c ^712) 



( V |yl e) |m)(c|^|n 2 )(n 1 n 2 |J(0)|t; / c) 



E c =Bl,-l V 711,712 

+ ^|vl e Vi)^ic|/(0)|^n 2 )(n 2 |^|c) + (^|W / |n 1 )(c|^ e) |n2)(nin 2 |J(0)K/c) 



(v| (mc|7(0)|^n 2 ) (n 2 |\4 e) |c) 



E 



711,712 



712/ L 



(c\Vi e) \ni){m\W\n 2 }{n 2 v\I{0)\cv') 



+ (c\Vg ) \n 1 )(n 1 v\I{0)\n 2 v')(n 2 \W\c) + (c|Vr|n 1 )(ri 1 |F] e) |n2>(^|/(0)|c ? ; / ) 



(SnT^Ec) 



E E 



(c| Vl e) I n> (n j W| ^ 1 7(0) | cu 7 ) 

£5=Bi,-i n 

+ (c\vi e) \n}(nv\I(0)\cv'}(c\W\c) + (c\vi e) \c}(c\W\n)(nv\I(0)\cv'} 



B 



(66) 
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r(3d) 



E 

ni,n2 



(^|vl e Vi)(c|Vr|r i2 )(n 1 n 2 |J(a;)| OT / ) 



(e w -£ m )(£ c -£„ 2 ) 



B 



+ (vl^lniXnicl/Hlnau'X^l^lc) 



+ E 

ni,ri2 



(^C ^ni) (^C ^712 7 

+ (c|vl e Vi)(^l^)l^ 2 )(n 2 |iy|c) + (clVTlniXrul^lnaXfiaul/Hlt/c) 



E E 



{c\Vi e) \n) {n\W\c) (cu\I{u) \v'c) 
+ (c\vi e) \n)(nv\I(uj)\v'd}(c\W\c} + (c\vi e) \c}(c\W\n) (nv\I(u)\v'c) 



B 



E E 



(c\vi e) \n)(n\W\c)(cv\I'(uj)\v'c) 
+ (c\vi e) \n)(nv\I'(uj)\v'c}(c\W\c) + (c\Vi e) \c){c\W\n)(nv\I'(oj)\v'c) 



E E 

ee=E\-\ n 

+ E E 

£{j — Si) Tl 

-E E 

£7} — &v £c — $1,- 



S v S 

2 



(v\vi e) \n)(nc\I'(u)\cv')(c\W\c) + (v\W\n){nc\I'(ou)\cu')(c\vi e) \c) 



(v\Vi e) \v)(vc\I'(cu)\nv')(n\W\c) + (v\W\v)(vc\l'(u)\nv')(n\vi e) \c) 



B 



<v|vl e) |C>{t)c|J /, (cc;)|5i/><a|Vr|c) + (^|PF|€;)(f;c|r / (o;)|^ , )(£| 1 lA ej |c) 



T"(,.M™A/Z\lA e )\ 



E E (c|^ e) |£)(^|rH|t/c)(c|^|c) 

££=^1,-1 £c=Sl,-l 



(67) 



Here 



1 for W = i4fs or H/ = 
i for H/ = 1/l e) 



r(Q) 



f 



(nin 2 \I(uj)\n 3 n4) = / df 1 df 2 w] 11 (x 1 )it] 12 (x 2 )J(a;)it n3 (fi)it re4 (^ 2 ) 

(1 - c?i • a 2 ) cos(a;ri 2 ) 



J(a;) = a- 



/'fa;) 



da; 



T\2 



d 2 I{uo) 
du 2 



(68) 

(69) 
(70) 
(71) 
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-» i a (2s) 

u = e v — E± t -i, and r\i = \X\ — x 2 \. As in the case of evaluation of <5(?hfs considered above, the 
diagrams corresponding to Eq. (16TT) are presented in Fig. 1 and the ones corresponding to Eq. 
( 1621) are presented in Figs. 3-5. Separating the interelectronic-interaction corrections is carried 
out according to Fig. [2] just as it was done in section [31 

The calculation of h±(k) was discussed in detail in Ref. [35]. We found that 



hi(k)B = Mp 



9j 



Y el (k) - sfjYV&k) 



VoB. (72) 



Calculating the other matrix elements, we obtain 

14 



h 2 (k)B 2 



(aZ) 



-U^\aZ)^ B) 2 /{m e c 2 



(73) 



u ' ; 2 1 + 1/2 



9j + 9 1 



-a*Z±{giS$\aZ) + {aZy^Q( 
h 2 {j, k) = 0. 



m e c\2 2/ + 3 {t) 

^2 ( aZ ) 



h 



21 



HqB 



Here the total functions are 



U$\aZ) = U 2 (aZ) + ^B s (aZ) + ^C s (aZ) + ... 
S$\aZ) = S 2 (aZ) + ^B^aZ) + ^C^(aZ) + . . . , 
T 2 {t \aZ) = T 2 (aZ) + ^B Q (aZ) + -^C Q (aZ) + ... 



(74) 
(75) 

(76) 
(77) 
(78) 



It must be stressed that the contributions of order 1/Z 2 and higher to these functions are not 
included in Eq. ( 1541) . Their evaluation would require consideration of some higher-order terms in 
operator (T4~H]) . The expansions flTUj) - (1751) are presented by anology with Eqs. (12"B"1) and (|2"7j) . The 
function 



U(hqB) 2 W 

determines the one-electron contribution. Calculations employing Eqs. ( 164 
interelectronic-interaction corrections Bg(aZ), B^aZ), and Bg(aZ) 



B s (aZ) 
B,(aZ) 

B Q (aZ) 



a 2 Z 3 m e c 2 
14(^) 2 
12 



(79) 

67l) yield for the 

(80) 



a 



2 ^ 9l Y^(F) ^BM F 
2160 



23a 4 Z 2 Q (^f) 2 Fi? c } (F) Mo5M F 



:82) 
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It must be noted that, because of the smallness of the contribution determined by Bq(o>Z), only 
Bg(aZ) and B^(aZ) were evaluated in the present paper. For checking purposes the calculation 
of these functions was performed in both Feynman and Coulomb gauges. The results of both 
calculations coincide with each other. 

Solving equation (154"1) . we finally obtain for Mp = —I + 1/2, / — 1/2, 



AE mag (x) = AE^pg 



AE (2s) 

ai(l + e{)M F x + e 2 SfV 

m e c z 



4M P 

ci(l + + c 2 (l + 5 2 + M 2 F S 3 )x 2 



± 2V 1 + 2/ + l 



(83) 



Here 



^2 

Si 



5 2 
5 3 



7^i 2 Fk^)+«(^+i)] 



(25) 



2? 1 



e (t), ^ / 7 , 2 23Q /m e c\ 2 
^2 \OiZ) — (aZ) 



nog'j \ n J i(2i - 1 



Tt\aZ) 



11 



Vf(aZ) 



(azy 

21+1 



(2s) 



-a 2 Z- 



12( 



-a 2 Z 



1 



6(ft + <7/) I 



g'A^aZ) - (aZf^Q ( 

^ t} (^) + (^) 2 ^g( 



m e c\ 2 4i" 2 + 4i" + 3, 



ft 



1(21 - 1) 



T 2 (t) (aZ) 



m e c\2 21 + 3, 



T«(aZ) 



2/ 



1 ^ 7 s 23 



* + q'i 



360 



H J 1(21-1' 



T 2 (t) (aZ). 



For F' = I + \ and M F = ±(I + ±), in contrast to Eq. Q, we have 



(84) 
(85) 



(86) 
(87) 



AE„ 



AE- 



(2s) 
HFS 



± di(l + 771)2; + ?7 2 



A£ 



(28) 

HFS ^2 



m P cr 



(89) 



where 



771 = a 2 Z" 



^72 = £2 



6 (ft - 2/^) 



^(aZ) + (aZf^Q {^ff (aZ) 
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U^(aZ) 



(90) 
(91) 



and the "— " and "+" signs correspond to Mp = —(I + \) and Mp = I + |, respectively. 

If / = 1/2, the electrical quadrupole interaction vanishes and one can easily obtain for Mp = 0: 



AE mag (x) 



AR 



(2s) 
HFS 



AE- 



e 2 - 



(2s) 

HFS _2 



x 2 ±-^l + c 2 (l + 5 2 )x 2 



(92) 



with 



ft 



6(ft + ft^) 
11 



(93) 



For / = 1/2, Mp = ±1, the effect is described by formula (1891) with 



'''=i^r 2zs?,(ttZ) - (94) 



5 Numerical results 

In Table 1, we present the numerical results for the functions U 2 (aZ), Bg(aZ), U^ 1 [otZ\ ^(aZ), 
B^aZ), S^\aZ), and T 2 {aZ) (only for the isotopes with I > 1/2) defined by Eqs. (1791). (150ft, 
f!76|) . (1291 . flHTl) . f!77|) . and fl30j) . respectively, for the 2s state. All the values are calculated for 
the extended nuclear charge distribution. The root- mean-square nuclear charge radii (r 2 ) 1 / 2 were 
taken from Ref. [44]. For those elements for which no accurate experimental radii were available 
we employed the empirical expression [45] 

(r 2 ) 1/2 = 0.836A 1/3 + 0.570(±0.05) fm , (95) 

where A is the nuclear mass expressed in a.m.u. The calculations were performed using the dual- 
kinetic-balance (DKB) basis set method [46] with the basis functions constructed from B-splines 
[47,48]. The uncertainties of U 2 (aZ), Bg(aZ), S 2 (aZ), B^aZ), and T 2 (aZ) were estimated 
by adding quadratically two errors, one obtained by varying (r 2 ) 1 / 2 within its uncertainty and 
the other obtained by changing the model of the nuclear-charge distribution from the Fermi to 
the homogeneously- charged- sphere model. The uncertainties of the total functions U 2 \aZ) and 
S 2 \aZ) due to uncalculated second- and higher-order terms were estimated as the first-order 
correction (~ Bg(aZ)/Z and ~ B fl (aZ)/Z , respectively) multiplied by the factor 2/Z. The 
uncertainty due to uncalculated first- and higher-order terms in Eq. ( 1781) was estimated in a 
similar way. 

In Table 2, we present the individual contributions to the 2s gj factor for some Li-like ions 
with / ^ in the range Z = 6 — 32. The Dirac point-nucleus value is obtained by Eq. (13). The 
interelectronic-interaction (Ag int ), QED (A^qed), nuclear-recoil (Ag^c), and nuclear-size (A#ns) 
corrections are obtained as described in Refs. [32,33]. The nuclear-size correction was evaluated 
for the homogeneously- charged- sphere model if Z = 6 — 16 and for the Fermi model if Z = 20 — 32 . 
The nuclear polarization contribution to the 2s gj factor of light Li-like ions can be neglected [17]. 
The gj factor values given in Table 2 are used for calculations of the coefficients in the Breit - 
Rabi formula. 

In Table 3, the numerical results for the coefficients in Eqs. ([2]), Qj, (Till . ( 1531 . ( 1891) . and 
( 1921) are listed for some Li-like isotopes in the interval Z = 6 — 32. Since in all the cases under 
consideration the absolute value of the recoil correction to the bound-nucleus gi factor is smaller 
than KT 11 [12], we actually have in Eq. ©: g\ = ^gi- 
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6 Discussion 



The energy separation between the ground-state HFS components (F = 1 — 1/2 and F 1 = I + 1/2) 
of a lithiumlike ion can be written as [49] 

AEU = -a(aZ) --—m e c 



x ^ 2s \aZ){l - 6^(1 - e&^+x™] 

+ L B M(aZ) + -LcW(aZ) + ...}, (96) 



449 

'aZ) 4 + ... (97) 



where 

A™(aZ) = 2 ' 2 , (1 + 1 + yf? 1 = 1 + ~(^) ■ 

V ; (1 + 7 ) 2 7 (4 7 2 - 1) 8 V ; 128 

is the one-electron relativistic factor, is the nuclear charge distribution correction, 

is the nuclear magnetization distribution correction (the Bohr - Weisskopf effect), a^ad * s ^he 
QED correction, B^ 2s \aZ) and C' 2s '(«Z) determine the interelectronic-interaction corrections to 
the hyperfine structure. Therefore, the dimensionless variable x = fi B/AE^p S is of order of 
6/io-B /[«(«Z) 3 ^m e c 2 ]. The intervals of B and Z, for which x ~ 1, are of special interest (in the 
original paper [36] the fields with 0.1 ^ x =Sj 3 were considered to be intermediate). 

For the magnetic fields with the magnitude B ~ 1 — 10 T, that are generally used in this kind 
of experiments, Li-like ions with Z = 6 — 32 meet the requirement x ~ 1. For this reason, only 
such ions are presented in Tables 1-3. 

For ions with Z < 32, the electric-quadrupole corrections to the coefficients a±, c±, c 2 , and d± 
are either equal to zero, if / = 1/2, as in the case of 13 C 3+ , or by a factor of 10~ 3 — 10~ 4 smaller 
than the magnetic-dipole ones. This is due to an additional factor (aZ) 2 in the electric-quadrupole 
contributions compared to the magnetic-dipole ones in the equations for the hyperfine-structure 
corrections to the Breit - Rabi formula coefficients and small values of Q for \ow-Z ions. 

As one can see from Table 3, the corrections e\, 5i, 62, S3, and rji for Li-like ions are several times 
smaller as compared to the corresponding ones for the Is state of the same H-like isotopes [21]. 
However, they provide more precise determinations of the coefficients in the Breit - Rabi formula. 

For B = 1 — 10 T, an estimate of the terms of the third and higher orders with respect to B 
in Eq. ( |54^|) indicates that the contributions from these terms are negligibly small as compared 
to both magnetic dipole and electrical quadrupole corrections. However, it is very important to 
take into account and 772-B 2 if Z = 6 — 32. This is due to the fact that these terms are 

comparable with the other corrections to the Breit - Rabi formula considered and the less Z is, 
the more appreciable the contributions from 62-B 2 and ^B 2 become. One can see that for Li-like 
ions these terms are 10 - 10 3 times bigger as compared to the case of the Is state of the same 
H-like isotopes [21]. In the second-order approximation (1541 with respest to B, formulas ([2]), ( fTTl) . 
( 1831) . and ( 1921) do not contain B to a power higher than two under the square root. This is due to 
the fact that h 2 (F) = h 2 (F'). 

The Breit - Rabi formula for the 2s state contains AE^p S , and the coefficients in the formula 
and the corrections to them calculated above include the value of fi/fiN- The uncertainties of the 
nuclear magnetic moments indicated in Table 3, as a rule, do not include errors due to unknown 
chemical shifts which, in some cases, can contribute on the level of a few tenths percents. Thus, 
carrying out the experiments on the Zeeman splitting with the aforesaid accuracy could provide 
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the most precise determination of both AE^p S and jj/jJN- The corrections to the Breit - Rabi 
formula evaluated in this paper will be important for this determination. 
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Table 1: The numerical results for the extended-charge- nucleus values of functions U% (aZ), 
S { 2 t] {aZ), and T 2 {t \aZ) (for the ions with / 7^ 1/2). The values of (r 2 ) 1 / 2 are taken from Ref. [44]. 



Ion 


13Q3+ 


17Q5+ 


21 Ne V + 


25 Mgy+ 


33gl3+ 


43 Ca iY+ 


53 Cr .l+ 


Z 


6 


8 


10 


12 


16 


20 


24 


(r 2 ) 1 / 2 , fm 


2.461 


2.695 


2.967 


3.028 


3.251 


3.493 


3.659 


U 2 {aZ) 


0.998574 


0.997464 


0.996038 


0.994295 


0.989858 


0.984153 


0.977179 


B s {aZ) 


2.47400 


2.47359 


2.47305 


2.47240 


2.47070 


2.46848 


2.46568 


U^(aZ) 


1.41(14) 


1.31(8) 


1.24(5) 


1.20(3) 


1.144(19) 


1.108(12) 


1.080(9) 


S 2 (aZ) 


1.00306 


1.00545 


1.00854 


1.01235 


1.02218 


1.03513 


1.05145 


B^aZ) 


-1.60040 


-1.60757 


-1.61684 


-1.62825 


-1.65769(1) 


-1.69639(1) 


-1.74505(2) 


S [ 2 t] {aZ) 


0.74(9) 


0.80(5) 


0.85(3) 


0.88(2) 


0.919(13) 


0.950(8) 


0.979(6) 


T 2 (aZ) 




1.00448(2) 


1.00710(2) 


1.01051(3) 


1.01927(4) 


1.03083(7) 


1.04531(10) 


T 2 W (aZ) 




1.0(3) 


1.0(2) 


1.01(17) 


1.02(13) 


1.03(10) 


1.05(9) 



Ion 


61 N p+ 




73 Ge 29+ 


Z 


28 


30 


32 


(r 2 ) 1 / 2 , fm 


3.822 


3.964 


4.063 


U 2 (aZ) 


0.968938 


0.964342 


0.959428 


B s {aZ) 


2.46226 


2.46030 


2.45817 


U^(aZ) 


1.057(6) 


1.046(5) 


1.036(5) 


S 2 (aZ) 


1.07146(1) 


1.08296(1) 


1.09555(2) 


B^aZ) 


-1.80458(3) 


-1.83874(3) 


-1.87609(4) 


S [ 2 \aZ) 


1.007(5) 


1.022(4) 


1.037(4) 


T 2 (aZ) 


1.06287(11) 


1.07277(14) 


1.08359(13) 


T 2 {t \aZ) 


1.06(8) 


1.07(7) 


1.08(7) 
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Table 2: The individual contributions to the ground-state gj factor of lithiumlike ions with 
nonzero nuclear spin and the nuclear charge in the range Z = 6 — 32. The values of (r 2 ) 1 / 2 are 
the same as in Table 1. 



Ion 


13Q3+ 


17Q5+ 


33gl3+ 


43 Ca iY+ 


9d 


1.999680300 


1.999431380 


1.997718193 


1.996426011 




0.000130758(19) 


0.00017666(3) 


0.00036124(9) 


0.00045445(14) 


A#QED 


0.002319417(6) 


0.002319549(12) 


0.00232070(6) 


0.00232171(10) 


Ao (e) 


0.000000009 


0.000000016 


0.000000045(1) 


0.000000057(2) 


A#NS 


0.0 


0.0 


0.000000005 


0.000000014 


9j 


2.00213048(2) 


2.00192760(3) 


2.00040018(11) 


1.99920224(17) 



Ion 


53Q.31+ 


73 Ge ^+ 


9b 


1.994838064 


1.990752307 




0.0005485(2) 


0.0007397(4) 


A#QED 


0.00232304(15) 


0.0023270(2) 




0.000000069(4) 


0.000000093(9) 


A#NS 


0.000000035 


0.000000160 


9j 


1.9977097(3) 


1.9938193(4) 
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Table 3: The numerical values of the coefficients in Eqs. ©, ©, flU), (PD, (PI), and (152]) for 
Li-like ions with Z = 6 — 32. The values of fi/fiN an d Q are taken from Refs. [50] and [51], 
respectively. 



Ion 


13Q3+ 


17Q5+ 


33gl3+ 


43 Ca iY+ 


/ 


1/2 


5/2 


3/2 


7/2 


/i//ijV 


0.7024118(14) 


-1.89379(9) 


0.6438212(14) 


-1.317643(7) 


Q, barn 




-0.02558(22) 


-0.0678(13) 


-0.0408(8) 






0.00041256(2) 


-0.0002337573(5) 


0.0002050317(11) 


ei 




-0.0000284(18) 


-0.0000653(9) 


-0.0000843(7) 


ai(l + ei) 




0.00041254(2) 


-0.0002337421(6) 


0.0002050144(11) 


e 2 {= V2) 


1.03(10)xl0 4 


5.4(3) xlO 3 


1.17(2)xl0 3 


7.28(8) xlO 2 


Cl 




2.00151505(4) 


2.00063394(11) 


1.99899721(17) 


*i 




0.0000000059(4) 


-0.00000000763(11) 


0.00000000864(7) 


Cl(l + 5i) 




2.00151506(4) 


2.00063393(11) 


1.99899722(17) 


c 2 


4.01159069(8) 


4.00606248(15) 


4.0025362(4) 


3.9959898(7) 


5 2 


-0.0000000151(18) 


0.0000000117(7) 


-0.0000000152(2) 


0.00000001730(15) 


5 3 




0.0 


-0.00000000001 


0.0 


c 2 (l + 5 2 ) 


4.01159062(8) 


4.00606253(15) 


4.0025361(4) 


3.9959899(7) 


c 2 5 3 




0.0 


-0.00000000002(1) 


0.0 


di 


1.000682697(10) 


1.00199519(5) 


0.99984946(5) 


1.00031873(9) 


m 


0.0000000075(9) 


-0.0000000292(18) 


0.0000000229(3) 


-0.0000000605(5) 




1.000682704(10) 


1.00199516(5) 


0.99984948(5) 


1.00031867(9) 



Ion 


53 Cr ^l+ 


™Ge 29+ 


/ 


3/2 


9/2 




-0.47454(3) 


-0.8794677(2) 


Q, barn 


-0.150(50) 


-0.196 


a\ 


0.000172295(11) 


0.00010643846(2) 


ei 


-0.0001041(6) 


-0.0001472(6) 


ai(l + ex) 


0.000172277(11) 


0.00010642279(7) 


e 2 (= V2) 


4.93(4) xlO 2 


2.660(13)xl0 2 


Cl 


1.9975374(3) 


1.9937129(4) 


Si 


0.00000000893(5) 


0.00000000776(3) 


Cl(l + 5l) 


1.9975374(3) 


1.9937129(4) 


c 2 


3.9901556(10) 


3.9748910(17) 




0.00000001803(11) 


0.00000001582(6) 


S3 


-0.00000000004(2) 


-0.00000000001 


c 2 (l + 62) 


3.9901557(10) 


3.9748911(17) 


C2S3 


-0.00000000018(6) 


-0.00000000005 


d\ 


0.99911329(13) 


0.9973886(2) 




-0.00000002699(17) 


-0.0000000708(3) 


d 1 (l + rj 1 ) 


0.99911326(13) 


0.9973886(2) 
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Figure 1: The second-order diagrams contributing to S%\aZ), T^(aZ) (if / = 1 and W 
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Figure 2: C is the original contour of the integration over the electron energy variable in the 
formalism with the standard vacuum. C is the integration contour for the vacuum with the (Is) 2 
shell included. The integral along the contour C- int = C — C describes the interaction of the valent 
electron with the (ls) 2 -shell electrons. 
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Figure 3: The third-order diagrams contributing to S^ ^aZ) (if / = 1 and W = Vjjpg) and 
U2 \aZ) (if / = ~ and = vl 6 '') being combined with products of the lower-order diagrams 
presented in Figs. 0] and EJ 



V 



(/0 



HFS 
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Figure 4: The first-order diagrams contributing to S^ictZ) and U^i&Z) being multiplied by the 
second-order diagrams presented in Fig. EJ 
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Figure 5: The second-order diagrams contributing to (qZ) and iptZ) being multiplied by 
the first-order diagrams presented in Fig. HI 
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